CASE FI 
COPY 


i 


X S2 64620 



92 1st 


• • • jCopy 

••• *•* RM E57K25a 



RESEARCH MEMORANDUM 


A METHOD FOR DETERMINING TURBINE DESIGN CHARACTERISTICS 

FOR ROCKET TURBODRIVE APPLICATIONS 

By Warner L. Stewart, David G. Evans, and Warren J. Whitney 

Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 









NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


WASHINGTON 

May 2, 1958 


CONFIDENTIAL 


• • 
• • 


NACA RM E57K25a 


••••• • ••• 

• •• • • • • • • 


CONFIDENTIAL 


• • • • • 

• • • 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


RESEARCH MEMORANDUM 


A METHOD FOR DETERMINING TURBINE DESIGN CHARACTERISTICS 

FOR ROCKET TURBODRIVE APPLICATIONS 

By Warner L. Stewart, David G. Evans, 
and Warren J. Whitney 


SUMMARY . 

This report presents a method for making a rapid determination of 
the design characteristics of turbines for rocket turbodrive applica- 
tions . Factors considered include relations between turbine flow and 
weight and rocket weight , as well as such turbine parameters as stage 
number, size, and pressure ratio. Independent variables used in the 
report include the amount of pump flow bled off for driving the turbine 
and the number of turbine stages . The design considerations presented 
herein are applied to an example rocket application using hydrogen and 
oxygen as the propellant to illustrate the type of results obtained 
through use of this method. 


INTRODUCTION 

The turbopump component of a rocket has the function of increasing 
the pressure level of the fuel and oxidizer from that of the tanks to the 
high pressure needed in the thrust chamber. The turbopump unit consists 
of the pumps which develop the pressure, the turbine that supplies the 
required power to the pumps, and the necessary gearing that connects the 
pumps and turbine. 

The ultimate merit of a turbopump system for a given rocket is in- 
dicated by minimum rocket gross weight. A consideration of rocket 
gross weight as affected by the turbopump system includes the effect of 
the dry weight of the components as well as the effect of the propellants 
consumed by the turbine. 

The optimization of the turbine component in this manner for any 
given application is complicated because a number of variables and l imi ts 
must be considered. Because of this it is important to have a design 
method available that includes these effects. The method should permit 
a rapid screening of a large number of turbine design’s to determine the 
configuration that best suits the application. 


CONFIDENTIAL 





2 


CONFIDENTIAL 


NACA RM E57K25a 


This report presents one such method. Independent variables consid- 
ered include the ratio of turbine flow to pump flow and the number of 
turbine stages. The report first considers equations to relate the tur- 
bine flow and weight variations to over-all rocket weight characteristics. 
Since the turbine weight is related to its frontal area in addition to its 
stage number, a subsequent part of the report is devoted to the develop- 
ment of equations necessary to obtain this quantity. In selecting the 
range of stage number, consideration is given to its dependency on rotor 
blade speed and specific work output. The analysis results of reference 
1 are used to do this. The necessary equations to obtain the turbine 
pressure ratio are included. The effect of turbine rotor disk and blade 
stress limitations on the selection of certain turbine parameters is also 
described. 

In order to illustrate the method and type of results that can be 
obtained, an example application will be described. This example case 
is a two-stage rocket using hydrogen and oxygen as the propellant. These 
fluids were chosen as representing an advanced application problem. Also, 
in the example, the weight of all components other than the turbine is 
assumed constant and therefore independent of variations in turbine char- 
acteristics. Under this assumption a range of turbine designs yielding 
minimum rocket gross weight is obtained. 


SYMBOLS 

A area, sq ft 

^ aspect ratio 

Cp specific heat at constant pressure, Btu/lb/°R 

e disk taper parameter 

F thrust, lb 

G ratio of thrust to gross weight 
g acceleration due to gravity, 32.17 ft/sec 2 
AH' total enthalpy requirement of turbine, Btu/sec 
h specific enthalpy, Btu/lb 

I specific impulse, sec 

J mechanical equivalent of heat, 778.2 ft-lb/Btu 
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K constant of proportionality in turbine weight equation (eq. (14)) 
n number of turbine stages 
P payload, lb 

p pressure, lb/sq ft 

R gas constant, ft/°R 

r radius, ft 

S structural parameter 

S structural parameter excluding effect of turbine weight 

s stress, lb/sq ft 

T temperature , °R 

t burning time, sec 

U turbine blade speed, ft/sec 

V absolute velocity, ft/sec 
W weight, lb 

w weight-flow rate, lb/sec 
y ratio of turbine flow to total pump flow 
y ratio of specific heats 

e ratio of specific impulse in exhausting turbine flow to that of the 
flow passing through the thrust chamber 

t) efficiency 


q g turbine efficiency based on static- to total-pressure ratio 
X work-speed parameter, U^/gJ Ah' 

p density, lb/cu ft 

$ rocket trajectory angle measured from vertical, deg 
oj rotative speed, radians/sec 
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Subscripts : 
an annulus 

B burnout 

b blade 

cr critical 

d disk 

e empty 

F frontal at turbine exit 

f fuel 

g gross 

H horizontal 

h hub 

id ideal 

2 last stage 

m mean 

o oxidant 

P pump or pumps 

p propellant 

ST stage 

s static 

T turbine 

t tip 

u tangential component 

v vertical 
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x axial component 

1 turbine inlet 

2 turbine exit 
Superscript: 

' absolute total state 


RELATION BETWEEN TURBINE AND ROCKET 

As pointed out in the INTRODUCTION, this report is concerned with a 
method for screening turbine designs for a specified turbopump applica- 
tion with the objective of obtaining a turbine design that is optimum 
from a rocket weight standpoint. Thus, it is necessary to first deter- 
mine equations that relate the two turbine parameters that affect rocket 
weight, turbine flow and weight, to the weight of the rocket. These equa- 
tions are now described. 


Relation Between Turbine Flow and Ratio of Rocket 
Gross Weight to Empty Weight 

Figure 1 presents a schematic diagram of a typical turbopump system 
where certain percentages of the fuel and oxidant are bled off at the 
pump exits. These fluids are then -burned in a gas generator, passed 
through the drive turbine, and thereafter are exhausted overboard. 

If none of the pump flow was bled off for powering the turbine, the 
thrust that would be obtained for a given thrust-chamber specific impulse 
I would be 


F = IWr 


( 1 ) 


However, since a certain percentage of the flow must be diverted through 
the turbine, this equation must be modified to 


= I _w, 


(1 - 6 ) 


v T 


or 


F = lw p [l - (1 - e) y J 


( 2 ) 
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where y = Wip/wp and represents the percentage of the total pump flow 

bled off at the pump discharge and e is the ratio of specific impulse 
achieved in exhausting the turbine flow overboard to that achieved in 
the thrust chamber. 

The manner in which y affects the rocket weight can he deter- 
mined through consideration of its trajectory and acceleration charac- 
teristics. Figure 2 presents a schematic diagram of a rocket in flight. 
Excluding the effects of aerodynamic drag and variation of acceleration 
of gravity with altitude, the following equations can he written relat- 
ing the force acting on the rocket and the acceleration along the verti- 
cal and horizontal directions: 


W V - L v \r 

F cos § - W = — 


dVy 

g dt 


y dVtT 
F sin $ = — • — i 
g dt 

Then, since 

W = Wg - wpt 

equations (2) and (4) can be substituted into equations (3) to yield 
Iw p [l - (l - e)y] COS § 


Wg - w p t 


1 _ 

g dt 


Iwp[l - (l - e)y]sin p dV H 
Wg. - wpt g dt~ 


(3a) 

(3b) 

( 4 ) 

(5a) 

(5b) 


In order to integrate equations '(5a) and (5b) I, e, and <t> must be 
known as a function of time. For simplification in this report, however, 
let I and e be assumed constant and two limiting trajectory paths be 
used. For the specifications of vertical flight, $ » 0°. Thus, equation 
(5a) can be integrated to 


AVi 


B 


gl|l - (1 - e)y] 


= in 



i[i' - (1 - e)y] 


(6a) 


If the trajectory is horizontal, # = 90° and equation (5b) can be inte- 
grated to 


AV, 


B 


gl[l - (1 - e)y] 


= In 


(6b) 
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Acceleration characteristics permit another equation to be developed 
which can be used in conjunction with the trajectory equations. By de- 
fining a parameter G as the ratio of thrust to gross weight, then 

Iw p [l - (1 - e)y] 

o«f 5 

w g w g 

Thus, 

Wp G 

Vi[i-U- 6 )y] 

Also, 


W g = W e + w p t B 


Dividing through by Wg gives 




Substituting equation ( 7 ) into (8) then yields 


( 8 ) 


W P 

1 = — + G 


W 


g 


H 

I[1 - (1 - e)y] 


or solving for Wg/W e gives 


W g 


% 


-1 


l[l- (l - e)y] 


(9) 


Equations (6) and (9) can be solved simultaneously for W g /W e and 
tp for given rocket requirements and limits over a range of the first 
independent variable y. Thus, the effect of varying y on gross-to- 
empty weight can be obtained. An example solution of equations (6) and 
(9) is presented in figure 3 in dimensionless form for G = 1.4 • The 
subsequent effect of variations in Wg/W e on the over-all rocket weight 
for a given payload is now described in conjunction with turbine weight 
considerations. 
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Relation Between Turbine Weight and Rocket Gross Weight 

The variations in rocket gross weight are studied in terms of the 
ratio of gross weight to payload W g /P using a structural parameter S. 
By definition, 

W p - P 
S = — 

*g - w e 


that is, the ratio of structure weight to propellant weight. Dividing 
through by W e gives 


S = 


1 . 

W g W e 
W e 

-£ - 1 

w„ 


Then solving for W g /P yields 



( 10 ) 


Figure 4 presents equation (10) in graphical form where W g /P is 
presented as a function of W g /w e for varying S where the effects of 
varying these parameters on gross weight can be determined. 


The derivations presented thus far are for one rocket stage. The 
effect of multistaging the rocket can be seen through consideration of a 
two-stage rocket as an example. In this case, the gross weight of the 
second stage is the payload of the first stage. Since primary interest 
is in the ratio of over-all rocket gross weight to second-stage payload 
(W g /p) over _ all , this parameter can be obtained from the equation 



over-all 



first stage 



second stage 


( 11 ) 


In the method of analysis presented herein, the effect of variations 
in turbine weight on the rocket weight is considered through variations 
in the structural parameter S. The parameter is comprised of the 
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contribution of all component weights. That is, 

S = (AS^ an ^. + ASp + ASg ear ^ n g + AS^^ers ) + ^®T 


= S + A£*ji (12) 

This report is concerned with, equations for obtaining Af>j. In an 

actual study, however, all the components must be considered to evolve 
an optimum compatible combination since, in general, the tanks, pumps, 
gearing, and turbine are interdependent. 

The equation for AStp is, by definition, 


w T w T 
^ = Wp 


(13) 


An equation for use in obtaining turbine weight was formulated from 
both published (ref. 2) and unpublished data and is expressed in terms 
of turbine frontal area representable by that at the turbine exit Ap 
and stage number n. The equation, as used herein, is 


Wrp 


Kh 1 / 2 A F 


(14) 


From single-stage tin-bine weight studies presented in reference 2, 
turbine weight was found to vary approximately as Ap for a value of K 
of approximately 70. Additional unpublished turbine weight studies veri- 
fied this value for K in addition to fixing the staging effect on tur- 
bine weight at approximately nl/2. 


Substituting equation (14) into (13) then yields 


AS 


T “ 


Kh 1 / 2 A F 
*B W P 


(15) 


Since tg is obtained from the rocket requirements and n is given as 
the second independent variable, the only other quantity needed to com- 
plete the calculation of ASip is the turbine frontal area parameter 
Ap/wp. The method used for obtaining Ap/wp is presented in the section 

DETERMINATION OF TURBINE FRONTAL AREA PARAMETER. The computed value of 
ASip can then be used to obtain S (eq. (12)). This quantity can finally 
be used together with Wg/W e obtained from the rocket trajectory and ac- 
celeration equations to compute the ratio of gross weight to payload 
(eqs. (10) and (ll)). 
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RELATION BETWEEN TURBINE AND PUMP 


Two turbine characteristics are related to the pump, namely power 
and rotative speed. The turbine power requirement can be obtained 
through consideration of the pump pressure requirement. 


The power required for the fuel pump is 

£p£ 


J AHj = w f 


Pf^Pjf 


and that for the oxidizer pump is 


Ap 1 

J ^ - w o 

o 'P,o 


The total power that the turbine must deliver is then the sum of these 
two powers; that is, 


J AH' 


APf 

Wf* 

p f iip,f 


+ M r 


p o\,o 


Dividing through by w p = Wp + w Q gives 

J AH' _ w f A Pf v p ^>o 

w p " w f + w Q PfHp^ f w f + w D P 0 Ti p ^ 0 

or 

^o 

ATT 1 _ 1 ^Pf w f A Pq 

Wp - Wo JPfTlp f + Wo JP 0 Hp 

M f W f 


(16) 


Thus, from equation (16) the total turbine power requirement can be ob- 
tained in terms of unit pump flow. This requirement is one of the quan- 
tities needed to determine the turbine characteristics. 


The other characteristic linking the turbine to the pumps is the 
rotative speed if no gearing is used. Variation in turbine rotative 
speed can be considered in terms of variations in the turbine frontal 
area parameter Ap/wp. Multiplying numerator and denominator of the pa- 
rameter by o> 2 yields 


Aj, «r|n 2 itU 2 

W P wpr 2 wpjo 2 


( 17 ) 
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Now, 


rpm = 3Cb/jt (l8) 

Solving equation (17) for cd and substituting into equation (18) give 

rpm = 30U t 

Equation (19) shows that for given turbine parameters (Ut and Ap/wp) 

an actual pump flow is required to obtain the turbine rotative speed. 

If direct drive is used between the turbine and either or both of the 
pumps, equation (19) must be used to determine Ap/wp as a function of 

U-t- If gearing is used, however, the value of Ap/wp can be selected 
primarily from turbine considerations. 



TUKBINE STAGE NUMBER CONSIDERATIONS 


As pointed out in the INTRODUCTION, turbine stage number is being 
considered an independent variable. The range of stage number as well 
as the associated efficiency is, however, dependent upon the particular 
value of turbine flow y selected. The manner in which this range of 
n as well as the associated required static- to total-pressure ratio 
across the turbine can be obtained is now presented. 


Reference 1 presents a method for obtaining the turbine efficiency 
based on the static- to total-pressure ratio for given turbine stage 
number, power, and speed requirements. Figure 5, a replot of figure 5 
of reference 1, shows the variation in this efficiency with a work-speed 
parameter X, defined as the ratio of the turbine mean section rotor 
blade speed squared to the required specific work output. That is, 


gj Ah' 


( 20 ) 


The parameter shown in the figure is the stage number n. If X is 
known, the associated range of n can be obtained. The upper limit on 
n (where limiting efficiency occurs) is obtained from the equation 

n = l/X (21) 

which represents a stage X of unity (see ref. l) . 
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Turbine Required Work- Speed Parameter 

The determination of the work- speed parameter X involves first 
the calculation of the turbine mean section blade speed U m and the 
specific work output Ah’ . In obtaining the equation for U m it is as- 
sumed that this speed is constant through the turbine and is equal to 
that at the exit (station 2). This mean section blade speed is related 
to the turbine-exit hub-tip radius ratio and tip speed by the equation 



As will be discussed in the section HUB-TIP RATIO AND TIP SPEED 
SELECTION, the values of (r^/r^) g and must be compatible with 

rotor disk and blade stress limitations. 

The turbine specific enthalpy required is 


Ah' 


AH 1 

w T 


or, if y is introduced. 


AH* 



y 


(23) 


The parameter AH’/wp is obtained from equation (16), and y is one of 

the two independent variables. Equations (22) and (23) can be substi- 
tuted into equation (20 ) to give 



(24) 


For the range of y selected, values of X can be obtained from 
equation (24) for use in determining the range of n from figure 5. 
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Turbine Required Pressure Ratio 


Figure 5 is used not only to obtain the variation in n but also 
to obtain the associated turbine static-to-total efficiency t] s .' The 
required turbine static- to total-pressure ratio can then be obtained 
using the following definition of ti s : 


^s 


Ah' 


Ah! . 
id,s 


(25) 


where Ah’ is the actual enthalpy drop (eq. (23)) and Ah! . is the 

ideal specific enthalpy drop corresponding to the static- to total- 
pressure ratio across the turbine. Since Ahj^g is related to this 

pressure ratio by the equation 


Ah 


id,s 


c T' 
P 1 



(r-i)/r 


(26) 


equations (23) and (26) can be substituted into equation (25) with the 
resultant equation solved for P 2 /p{ to yiel^ 



\r/(r-i) 

AH* \ 

V P 

y c P T iv 


(27) 


It might be noted that, if, for a given application, a lower limit 
on P 2 /pi is specified (see example), equation (27) can be used to de- 
termine the resultant variation in t] s with y (and, hence, X) . in 

this manner a lower limit on n can also be obtained when these limiting 
calculation results are drawn on figure 5. 


DETERMINATION OF TURBINE FRONTAL AREA PARAMETER 

As pointed out in the section entitled Relation Between Turbine 
Weight and Rocket Gross Weight, the turbine frontal area parameter Ap/wp 

must be known or obtained before the rocket weight calculations can pro- 
ceed. If the turbine is given flexibility of speed, the determination 
of this parameter is made through consideration of the required turbine 
pressure ratio and continuity at the turbine exit (station 2). By 
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assuming uniform flow conditions, this continuity equation can be written 
as 


or 


W T - ( pV x) 2^&n,2 

= (p'V cr ) A an,2( p ' V cr)2 


= ( 7M a (P ’ Vcr) 2 ( p * V ) 

Ip'j U^j an >2 (p'V ) cr 1 

v / 2 ' cr '7> V cr'i 




P* 



w = /jlV^YM a/i ei § 

T \^>7 2 \v cr y V r + 1 R 



(28) 


Now, 


Also, 


(T(C • (« 1 ■ « 


- ■ (H 


'an 


(29) 


“ Att 


& 


( 30 ) 
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In addition, 

Ah ' = c p ( T i - 


A 




Solving for T^/T^ and introducing equation (23) yield 



Substituting equations (29), (30), and (3l) into equation (28) gives 



Dividing both sides by Wp and solving for Ap/wp result in 



One quantity in equation (32) requiring further derivation in this sec- 
tion is the temperature ratio (T/T’Jg* From conservation of energy. 
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( 33 ) 


As will be 
(Yx/^cr^2 

tangential 


discussed, the axial component of critical velocity ratio 
is selected through limiting-loading considerations . The 

component (V u /V cr )2 is considered as follows: 



Vu,2 AV u,I U m 
U m Y CTf2 


(34) 


By using equation (22) and the equation for ^cr-> an equation for the 
last term in equation (34) U m /V cr g can be written as 



or, by using equation (3l), 



and can be thus obtained. 

The parameter AV u ^/u m is by definition equal to the reciprocal 

of the stage work-speed parameter. Thus, using equal stage work as was 
done in reference 1 gives 


u 


m 



(36) 
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Finally, the parameter V Uj 2 /aVu, 2 is obtained from the blade velocity 
diagram considerations . From reference 1 the following conditions are 
assumed to prevail: 

For 0 < \grp < 0.5, 


V u,2 1 

AV u,2 " 2 



For 0.5 < < 1* 


u, 2 


AV, 


= 0 


u,2 


(37a) 


(37b) 


Thus, by using equations (34) to (37) the temperature ratio (t/T') 2 can 
be obtained by equation (33). . 

The selection of a value of (V x /V cr ) 2 must be made with consider- 
ation of the turbine limiting- loading characteristics. If a turbine is 
operated at a specified speed, decreases in static- to total-pressure 
ratio result in increases in specific work output. This occurs until a 
certain pressure ratio is reached where maximum work output is obtained. 
Further decreases in pressure ratio then reduce the efficiency without 
increasing the work output. 


The pressure ratio where maximum work is first obtained is termed 
the limiting-loading point and is a result of choking in the plane of 
the rotor trailing edge at the turbine exit. With no trailing-edge 
blockage the point would occur at the same time as choking with the exit 
annulus at station 2. However, as a result of this blockage in the ac- 
tual case, this point is obtained at values of exit axial Mach number on 
the order of 0.7. Since the usual practice is to design the turbine with 
some margin with respect to limiting-loading, axial Mach numbers less 
than 0.7 are usually prescribed. These considerations must then be made 
in selecting a value of (V x /v cr ) 2 . 


HUB-TIP RATIO AND TIP SPEED SELECTION 

From equations (24), (32), and (35) it is evident that values of 
hub-to-tip radius ratio and tip speed at station 2, (r^/ r^.) 2 and U^. 2 , 
must be selected before the aerodynamic aspects of the turbine design 
study can be considered. The selection of these two parameters must be 
made with consideration to disk taper and rotor blade hub centrifugal 
stress limitations. These considerations will now be made. 
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Turbine Disk Limitations 

The turbine disk limitations are considered herein through consider- 
ation of the taper of the disk from the rim to the disk center. This 
taper is illustrated in figure 6. Reference 3 describes this problem in 
detail and develops an equation relating the disk taper to the turbine 
hub-tip radius ratio. This equation, which was obtained for constant 
disk stress and no center hole, can be written in terms of the symbols 
used in this report as 



Figure 6 presents equation (38) in graphical form. The disk taper 
parameter e is shown on the left as a function of the hub- tip radius 
ratio for nontapered blades and stress -density parameters p^s^/p^s^ of 

1 and 2/3. The plot is for an aspect ratio of 2. These curves can be 
used to determine whether or not the selected hub-tip radius ratio at the 
turbine exit lies within prescribed disk taper and stress limitations. 


Turbine Rotor Blade Stress Considerations 


Reference 3 also presents the development of an equation relating 
the turbine centrifugal stress at the rotor blade hub to the blade tip 
speed and hub-tip radius ratio. This equation (eq. (8) in ref. 2) can 
be rewritten in terms of this report's symbols for the last turbine 
stage as 



where s-^ is the rotor hub centrifugal stress using a nontapered blade. 
Tapering the blade would reduce the stress somewhat. 
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Equation (39) is also presented in graphical form in figure 6. The 
tip speed parameter ( is shown at the right as a function of the 
hub-tip radius ratio r-^/r^. . From this curve it can be determined if the 
selected values of tip speed and hub-tup radius ratio at the turbine exit 
are within the rotor blade stress limitations. 


APPLICATION TO EXAMPLE ROCKET 

To illustrate the type of results obtained when applying the equa- 
tions described in this report, an example application will be described. 
It will be assumed for this application that variations in the turbine 
characteristics do not affect the weight of other components, thereby 
permitting a constant value of structural parameter, excluding the tur- 
bine, to be used. Furthermore, the turbine will be given the flexibility 
of rotative speed and will use a tip speed and hub-tip radius ratio 
determined from limiting rotor disk and blade characteristics. Under 
these assumptions and conditions the rocket gross weight characteristic 
will be presented over a range of y and A p /wp to indicate the area 
of design that yields the lowest rocket gross weight. 

The rocket will be assumed to have the following specifications: 

(1) Rocket requirements and assumptions 

(a) Vertical flight 

(b) Two stages 

(c) Final burnout velocity, 25,000 ft/sec 

(d) Change in velocity per stage AVg ^ ST , 12,500 ft/sec 

(e) Thrust to gross weight, G = 1.4 

(f) Zero aerodynamic drag 

(2) Structural parameter without turbine weight, S = 0.08 (both 
stages) 

(3) Propellant considerations (both stages) 

(a) Hydrogen fuel (p p = 4.42 lb/cu ft) 

'(b) Oxygen oxidant (p Q = 71.2 lb/cu ft) 
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(c) w o /w f = 3.2 

(d) Specific impulse, I = 360 sec 

(4) Pump 

(a) App = 700 lb/sq in. = 100,800 lb/sq ft 

(b) T]p = 0.60 

(5) Turbine (using fuel rich mixture) 

(a) U t ^ 2 = 1400 ft/sec 

(b) (r h /r t ) 2 = 0.80 

From figure 6 these quantities are compatible with 
s-jj = 40,000 psi 
= 526 lb/cu ft 
e = 4 

Pa s b/pb s a = 1 

(c) c = 2.12 

p 

(a) r = i.36 

(e) R = 437 

(f) T^ = 1860° R 

(g) = 0.90 Ap^ = 90,720 lb/sq ft 

(h) (V x /V cr ) 2 = 0.40 


(i) Turbine weight constant, K = 70 


(j) (p_/p') . , = 0.05 

2 1 minimum 


(k) £ = 0 
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Determination of Turbine Frontal Area Parameter 

Before calculation of the rocket weight characteristics can be 
made, variation in the turbine frontal area parameter Ap/wp over the 

range of selected flow y and stage number n must be determined. The 
procedure used to obtain this is outlined in table I together with typi- 
cal calculation results obtained for the example case. Columns 1 to 9 
are the quantities either specified or computed by the indicated equa- 
tions and represent input constants for the turbine aerodynamic phase of 
the study. 

Column 10 represents the first independent variable being consid- 
ered, y. Column 11 is the computed value of X required before enter- 
ing figure 5 to obtain the range of stage number n. Column 12 is then 
the other independent variable n, with the. upper limit being l/X and 
the lower li m it occurring at the limiting turbine pressure ratio. Fig- 
ure 5 is also used to obtain the turbine required static-to-total effi- 
ciency t] s shown in column 13. Columns 14 and 15 finally represent the 
calculations necessary to obtain the turbine frontal area parameter 
A F /wp (column 15). The calculations necessary to obtain each column are 
indicated at the bottom of table I where the quantities needed and equa- 
tions used are specified. 

The actual plot used to obtain n and tj is presented in figure 
7. The circles shown in this figure represent the points at which the 
calculations were made. 


Determination of Rocket Weight 

Once the turbine frontal area parameter is known, the rocket weight 
calculations can commence. Columns 16 to 21 of table I itemize the addi- 
tional required rocket and turbine constants. Then columns 22 to 26 pre- 
sent the necessary calculations to obtain the ratio of over-all rocket 
gross weight to payload (column 26). 

The results of the calculations made for the example case are pre- 
sented in figure 8. This figure shows a composite of the rocket and tur- 
bine primary parameters over the range of turbine flow y considered. 

The turbine parameters include the frontal area parameter Ap/wp, pre- 
sented as the* ordinate , with the stage number n shown as solid lines 
and the static- to total-pressure ratio shown as dot-dashed lines. The 
upper crosshatched area represents the specified lower limit on this 
pressure ratio. The lower crosshatched area represents the limitation 
imposed by the attainable efficiency (fig. 7). Any point within this 
bounded area represents a turbine design that will meet the pump power 
requirements . 
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Also shown in this figure are contours of computed over-all rocket 
gross weight to payload (W g /p) over _ all . Contours of 70, 68, and 67 as 

well as a line indicating the minimum value of 66.5 are included. * In- 
spection of 'this figure shows that, for this example case, the turbine 
flow has the predominate effect on the change in ( Wg /P ) 0 vkr-all 

the effect of Ap/wp being of secondary importance. This occurs be- 
cause, for a given y, , as Ap/Wp is varied the stage number n also 
varies in such a manner as to yield relatively small changes in turbine 
weight. It can also be seen that when consideration is given to the 
rocket weight characteristics a fairly well defined turbine design 
range emerges. One such turbine design could have the following 
characteristics : 

y = 0.01 

A F /wp = 17*icr 4 

n = 4 stages 

P 2 /p{ 15 0*05 

and from figure 7 

Ti =0.65 

■s 


CONCLUDING REMARKS 

This report has described a method of screening turbine designs for 
a rocket -pump drive application. In the example case mini mum rocket 
gross weight was used as a criterion, with all rocket components except 
the turbine assumed to be fixed or constant and therefore not affected 
by turbine operational characteristics. Rocket gross weight was used in 
this manne r to illustrate the combined effects of turbine weight and tur- 
bine propellant consumption on over-all rocket performance. In order to 
optimize an entire rocket design for a given application, the weights of 
the other major components such as propellant tanks, propellant pumps, 
and turbopump gearing (if any) must be considered, as well as their inter- 
dependent variation with one another. These additional factors must then 
be considered in an actual application to achieve an optimum compatible 
combination, which involves a study beyond the scope of this report. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, January 30, 1958 
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Figure 2. - Schematic diagram of rocket showing relation between thrust, 
weight, and attitude. 


CONFIDENTIAL 



[^(a ~ I) ~ ill pug 


[*(3 - T) - l!l8 


AV 


‘ s ja^sraejtBd auixx pue jCq.tooxaA 


_ CCSSTFIIiENTIAL 


Figure 3. - Rocket velocity- weight characteristics. Zero aerodynamic drag, no variation in g with altitude. 
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Ratio of gross weight to empty weight, W /w 

g e 

Figure 4. - Effect of structural parameter and ratio of gross weight to 
empty weight on ratio of gross weight to payload. 
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Figure 5. - Multistage turbine efficiency characteristics (ref. l). 
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Figure 6. - Turbine disk taper and rotor stress considerations. 
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Figure 7. - Application of multistage turbine efficiency characteristics to example case. 



Turbine frontal area parameter, Ap/vp, sq ft/(lb/i 
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Figure 8. - Example turbine characteristics including effect on ratio of over-all 
rocket gross weight to payload. 
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